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Connections of activated hopping processes with the breakdown of the Stokes-Einstein 
relation and with aspects of dynamical heterogeneities 

Song-Ho Chong 

Institute for Molecular Science, Okazakt 444-S585, Japan 

(Dated: July 27, 2008) 

We develop a new extended version of the mode-coupling theory (MCT) for glass transition, which 
incorporates activated hopping processes via the dynamical theory originally formulated to describe 
diffusion-jump processes in crystals. The dynamical-theory approach adapted here to glass-forming 
liquids treats hopping as arising from vibrational fluctuations in quasi-arrested state where particles 
are trapped inside their cages, and the hopping rate is formulated in terms of the Debye- Waller 
factors characterizing the structure of the quasi-arrested state. The resulting expression for the 
hopping rate takes an activated form, and the barrier height for the hopping is "self-generated" in 
the sense that it is present only in those states where the dynamics exhibits a well defined plateau. 
It is discussed how such a hopping rate can be incorporated into MCT so that the sharp nonergodic 
transition predicted by the idealized version of the theory is replaced by a rapid but smooth crossover. 
We then show that the developed theory accounts for the breakdown of the Stokes-Einstein relation 
observed in a variety of fragile glass formers. It is also demonstrated that characteristic features of 
dynamical heterogeneities revealed by recent computer simulations are reproduced by the theory. 
More specifically, a substantial increase of the non-Gaussian parameter, double-peak structure in 
the probability distribution of particle displacements, and the presence of a growing dynamic length 
scale are predicted by the extended MCT developed here, which the idealized version of the theory 
failed to reproduce. These results of the theory are demonstrated for a model of the Lennard- Jones 
system, and are compared with related computer-simulation results and experimental data. 

PACS numbers: 64.70.pm, 61.20.Lc 



I. INTRODUCTION 



The decoupling of the self-diffusion constant from the 
viscosity or the structural relaxation time - also referred 
to as the breakdown of the Stokes-Einstein (SE) rela- 
tion - that occurs for temperatures T < 1-2 Tg near the 
glass transition temperature Tg is among the most promi- 
nent features of fragile glass formers [l|, 0, Q. The de- 
coupling has been considered as one of the signatures of 
spatially heterogeneous dynamics or "dynamical hetero- 
geneities" [1, la, B 0] ■ O'^ the other hand, it is also well 
recognized that the onset temperature « 1-2 Tg of the 
decoupling is close to a crossover temperature at which 
transport properties change their characters |8i] and be- 
low which the dynamics is thought to be dominated by 
activated hopping processes over barriers [9|. Then, a 
natural question arises as to possible connections among 
the decoupling, dynamical heterogeneities, and the hop- 
ping processes. In this paper, such connections are ex- 
plored by extending the idealized mode-coupling theory 
(MCT) for glass transition [13]. 

The idealized MCT has been known as the most 
successful microscopic theory for glass transition. In- 
deed, extensive tests of the theoretical predictions car- 
ried out so far against experimental data and computer- 
simulation results suggest that the theory deals prop- 
erly with some essential features of glass-forming liq- 
uids [lliliJI- On the other hand, a well-recognized limi- 
tation of the idealized MCT is the predicted divergence of 
the a-relaxation time at a critical temperature Tc - also 
referred to as the nonergodic transition - which is not 



observed in experiments and computer simulations. An 
extended version of MCT developed in Ref. [13,] aims at 
incorporating activated hopping processes which smear 
out the sharp nonergodic transition and restore ergodic- 
ity for T < Tc, but its applicability has been restricted 
to schematic models. This is because of the presence of 
the subtraction term in the expression for the hopping 
kernel, which violates the positiveness - a fundamental 
property - of any correlation spectrum. 

There have been relatively few other attempts to go be- 
yond the idealized MCT [ij, [T^ [l3| , and incorporating 
the hopping processes into the theory for glass transition 
has been a major unsolved problem. We present here a 
new formulation which is motivated by ideas from the dy- 
namical theory originally developed to describe diffusion- 
jump processes in crystals [ITJ ]. The dynamical-theory 
approach adapted in this work to glass-forming liquids 
treats hopping as arising from vibrational fluctuations in 
quasi- arrested state where particles are trapped inside 
their cages, and the hopping rate is formulated in terms 
of the Debye- Waller factors characterizing the structure 
of the quasi-arrested state. The resulting expression for 
the hopping rate takes an activated form, and the bar- 
rier height for the hopping is "self-generated" in the sense 
that it is present only in those states where the dynam- 
ics exhibits a well defined plateau. It will be discussed 
how such a hopping rate can be incorporated to develop 
a new extended version of MCT. 

We will then investigate whether the developed the- 
ory accounts for the mentioned decoupling for T "^ T^. 
Such an investigation makes sense since Tc is also found 
to be close to 1.2 Tg [^, i.e., the onset temperature of 



the decoupling [l|, ^^ . It will also be examined whether 
our theory reproduces characteristic features of dynam- 
ical heterogeneities revealed by recent computer simula- 
tions. More specifically, we shall study whether the the- 
ory predicts a substantial increase of the non-Gaussian 
parameter [l9|, double-peak structure in the probability 
distribution of particle displacements [20, [211. and the 
presence of a growing dynamic length scale [6|, which 
the idealized MCT failed to reproduce. 

The paper is organized as follows. In Sec. |TT1 we for- 
mulate our new extended MCT. Numerical results of 
the theory are presented in Sec. IIIII for a model of the 
Lennard-Jones system, and connections of the hopping 
processes with the breakdown of the SE relation and with 
aspects of dynamical heterogeneities are discussed. The 
paper is summarized in Sec. IIVI Appendix [X] outlines 
a novel derivation of the hopping kernel formulated in 
Ref. [1^, and Appendix [B] is devoted to the derivation of 
the extended-MCT equations for the mean-squared dis- 
placement and the non-Gaussian parameter. 



II. THEORY 

We start from surveying basic features of the ideal- 
ized MCT 01 (see also Appendix E|. A system of N 
atoms of mass M distributed with density p shall be 
considered. Structural changes as a function of time t 
are characterized by coherent density correlators 4'q {t) — 
{pte''^*Pf^)/NSg. Here p^ = I]iexp(ig • n) with fi re- 
ferring to ith particle's position denotes density fluctua- 
tions for wave vector q; C the Liouville operator; (•) the 
canonical averaging for temperature T; Sq = {PffPq)/N 
the static structure factor; and q = \q\. Within the 
Zwanzig-Mori formalism |22l | one obtains the following 
exact equation of motion: 



Here p = q — k, and the vertex function is given by 



dUq{t)+nl<i>q{t) 



ni 



dt'mq{t-t')dt'(l)q{t') =0. 



(la) 

Here il^ = q^k-QT/MSq with Boltzmann's constant fce, 
and the memory kernel mq(t) describes correlations of 
fluctuating forces. Introducing the Laplace transform 
with the convention f{z) ~ iL dte'''^*f{t) (Imz > 0), 
Eq. (jlap is equivalent to the representation 



M^) = -i/{z-nl/[z + nlmqiz)]}. 



(lb) 



Under the mode-coupling approximation, the fluctuat- 
ing forces are approximated by their projections onto the 
subspace spanned by pair-density modes p^Pp- The fac- 
torization approximation for dynamics of the pair-density 
modes yields the following idealized-MCT expression for 
the memory kernel to be denoted as ■m^^{t): 



I' 



V{q;k,p) 



2{2Tifq- 



-SqSkSp[[q-k)ck + {q-p)CpY, (2b) 



in terms of Sq and the direct correlation function Cq = 
(1 - l/Sq)/p. The ideahzed-MCT equations (P) and (HD 
exhibit a bifurcation for (/'g(t — > oo) — fq - also referred 
to as the nonergodic transition - at a critical temperature 
Tc [l3|. For T > Tc, the correlator relaxes towards fq — 
as expected for ergodic liquid states. On the other hand, 
density fluctuations for T < T^ arrest in a disordered 
solid, quantifled by a Debye- Waller factor fq > 0. 

It is the factorization approximation [cf Eq. (jA7p ] that 
leads to the nonergodic transition at Tc. Therefore, one 
has to consider corrections, mq{z) — m^^ {z) + Amq{z) , to 
go beyond the idealized MCT, which shall be quantified 
via the hopping kernel defined by 



5q{z) - -l/mq{z) + l/mf{z) 



(3a) 



The correction term reads Amq{z) = m'T{z)dq{z)mq{z), 
and the memory kernel mq{z) can be expressed as 



mq{z) = n^^{z)l[l-5q{z)mf{z)] 



(3b) 



rr^^ it) = I dk V{q- k,p) <j>k (t) </>p (i) . 



(2a) 



As demonstrated in Appendix \^ one can derive based 
on Eq. (|3ap an expression for 5q{z) which is essentially 
the same as that of the extended MCT of Gotze and 
Sjogren [ij, [23| by applying the Zwanzig-Mori formal- 
ism to mq(t) and then introducing the corresponding 
mode-coupling approximation for the "memory kernel" 
to vriqit). Such a hopping kernel, however, retains the 
same problem mentioned in Sec. HI Instead, our at- 
tempt for the extension of the idealized MCT is moti- 
vated by the following observation: substituting Eq. (j3bp 
into Eq. ([Tb| yields for small z [11| 

^q{z) = -ll{z + Sq{z) -Ql/[z + nl mjf (z)]}. (4) 

Dropping Sq{z), this equation reduces to the one of the 
idealized MCT: approaching Tc from above, mjf (z) for 
small z becomes larger, and so does (t>q{z), leading to the 
nonergodic transition at T = Tc. In the presence of Sq{z), 
on the other hand, the transition is cutoff since the third 
term in the denominator of Eq. ^ becomes unimportant 
when m|f'(z) becomes large. The long-time dynamics of 
(j}q{t) in this case is thus determined by Sq{z) for small 
z. This observation raises a possibility of constructing a 
new approximate theory for Sq{z) described below. 

We first derive a rate formula for a hopping process 
in which an atom at r^ jumps to a nearby site r^ sep- 
arated by an interparticle distance. The presence of 
such a process at low temperatures has been revealed 
by computer simulations 24, 25, 26]. This will be done 
via the dynamical theory originally developed to describe 
diffusion-jump processes in crystals [17|. The approach 
adapted here to glass-forming liquids treats hopping as 
arising from vibrational fluctuations (phonons) in quasi- 
arrested state where particles are trapped inside their 



cages. Let us suppose that the quasi-arrested state char- 
acterized by the Debye-WaUer factors fq can be described 
as a frozen, irregular lattice [23 • Each particle then has a 
well defined equilibrium position Ri within the lifetime of 
the quasi-arrested state, and we introduce the displace- 
ment from the equilibrium position via fi — Ri + u{Ri). 
The essential feature of the hopping process is that a 
jumping atom passes over a barrier formed by neighbors 
which block a direct passage to the new site. The crite- 
rion that determines whether or not a given fluctuation 
is sufficient to cause a jump is therefore concerned with 
the relative displacements of the atom and the saddle- 
point. We thus employ as a "reaction coordinate" [rn |. 
x{t) = [u{Ri + s,t) — u{Ri, t)] ■ s, and assume that a hop- 
ping occurs when x{t) exceeds a critical value x*, which 
measures the size of fluctuation needed to cause a jump. 
Here s denotes the saddlepoint position with respect to 
Ri , and the scalar product selects only those fluctuations 
directed towards s = s/s. Each phonon displaces a hop- 
ping atom towards the saddlepoint. The phonon phases 
are random, but the displacements may occasionally co- 
incide in such a way that a hopping process occurs. The 
hopping rate Whop can then be calculated from such a 
probability per unit time, and one obtains along the line 
described in Ref. [I4I with the isotropic Debye approxi- 
mation, Whop = (l/27r)(3/5)i/2^Dexp[-3mi;2AV2fcBT], 
in terms of the sound velocity v. Here lod = k^v with 
the Debye wave number k^ — (67r^p)^/'^, and A = x* /s. 
Notice that the sound velocity here refers to the one in 
the quasi-arrested state, which is renormalized by the 
Debye-WaUer factors [28|, [2^. To emphasize this, the 
sound velocity shall be expressed as u = ^J Mi^/(pm) in 
terms of the (longitudinal) elastic modulus 



A/L = M{? + efL, 



(5a) 



consisting of the equilibrium value M^ — p{k'B,T)/ Sq, 
where Sq = 5g_,o i ^^nd an additional contribution for the 
quasi- arrested state, for which MCT yields [30| 



5Mi^^p{kBT) / dkVL{k)fl 



with 



^L(fc) 



47r2 



'dkc'klck 



[fee'" 



The hopping rate is then given by 



cjd cxp 



"^^-=2^15 



3Ml 
'2pkBT 



A^ 



(5b) 



(5c) 



(6) 



The pre-exponential factor represents a mean attack fre- 
quency while the exponential term of the activated form 
gives the probability that the system is found at the crit- 
ical displacement x*. In addition, the barrier height for 
the hopping is "self-generated" in the sense that it is de- 
termined by the plateau heights fq of the coherent density 



correlators, and is present only in those states where the 
dynamics exhibits a well defined plateau. 

We next relate the hopping rate Whop to the hop- 
ping kernel Sq(z). Our discussion becomes simpler if the 
tagged-particle density correlator (j)q{t) - the self part 
of (j)q{t) - is considered, so this case shall be consid- 
ered first. Hereafter, quantities referring to the tagged 
particle shall be marked with the superscript or sub- 
script "s". The Zwanzig-Mori equation for (j>q{t) has 
the same form as Eq. (fla| with 0^, Wg, and fl^ re- 
placed by (j)q, m^, and (flq)^ — q^k^T/M, respectively; 
the idealized-MCT kernel corresponding to Eq. ([2a|) is 
given by m^''i(t) = J dkV%q;k,p}4>kit)4>pit) with V = 
pS'fe[(g-fc)cfc]V[(27r)3g4] |3ij. and Eq. dSbJ holds with : 



,id 



'qj 



TO'", and Sq replaced by to^, m*'*^, and S^, respectively. 
The arrested part /| of the correlator 0^ (t) is referred to 
as the Lamb-Mossbauer factor. 

In the absence of the hopping kernel, the idealized ker- 
nel rriq "^ (t) for T < Tc arrests at a plateau for long times 
whose height is given by C| = /q /(I — fq) [lOl, i-e., there 
holds TO^"^(z) = ~Cq/z for small z. Substituting this 
into Eq. ([4]) for 4'q{z) yields for small z 

cl^'q{z)^-f^/[z + f^S^q{z)], (7) 

which determines the a relaxation of 4>g{t) - the decay 
from the plateau fq to zero - in the presence of Sq{z). 

On the other hand, when the a relaxation is donunated 
by hopping processes characterized by a rate Whop i'r — > 
r'), the van Hove self correlation function Gs{r,t), re- 
lated to (f>q{t) via the inverse Fourier transform 



G.(f,t) = 



(27r) 



dqe-^^^-UUt)> 



(8) 



and proportional to the probability of finding the tagged 
particle at r and t [22| , obeys a simple rate equation 

dtG,{f,t) - ^[wi,op(r + r^rOG,(f+M) 



-Whop{r^r + l)Gs{r,t)]. (9) 

We assume that only hoppings with \l\ ^ a are relevant, 
in which a = /q""° drr [iV(r)/iVc] denotes the weighted 
average of interparticle distances. Here rmin denotes the 
first minimum of the radial distribution function g(r) 
defining the first shell; N{r)dr with N{r) = Anr^ pg{r) 
gives the mean number of particles at distance between 
r and r + dr; and Nc = /n"'" drN{r) is the coordination 
number of the first shell. The quantity a serves as an 
analogue of the lattice spacing in crystals. Since there 
is no site dependence in the hopping rate we formulated 
[c/. Eq. ®], there holds 

9tG,(f,i) = u.hop 5I[G.(^+M)-G,(r,t)]. (10) 



Fourier transforming this yields 
9f^g(i) = -"'hop J2^^ 



-iq-l 



]m)- (11) 



Assuming that I are oriented at random, the summation 
^l,-, is given by the orientational average multiplied 

by the number of sites satisfying \l\ « a, which is ap- 
proximated by the coordination number N^ of the first 
shell. This leads to 

5t0^(t) = -WhopA^c[l - sin{qa)/{qa)](j>l{t). (12) 

Noticing that the a relaxation of 4)g{t) starts from the 
plateau /^ , the Laplace transform of this equation reads 

^li^) = -/,/{ 2 + *u.hopiVc[l - smiqa)/{qa)] }. (13) 

By comparing Eqs. ([7]) and (fT3| . one arrives at the fol- 
lowing expression for the hopping kernel: 



6^(z) 



iwhopNc[l - sin{qa)/{qa)]/f^ 



(14) 



The collective hopping kernel Sq{z) consists of the self 
and distinct parts, Sq{z) = Sg{z)/Sq + (5q"^*(z). In the 
present work, we shall adopt a simple model for 6q{z) 
in which the distinct part describing possible correlated 
jumps is neglected: 



Sq{z)^S^qiz)/Sq. 



(15) 



This model for 6q{z) looks oversimplified, but nontrivial 
theoretical predictions follow from such a simple model 
as will be demonstrated in Sec. IIIII 

Equations ^ and ^ with Eqs. (gal, §^, (fU)) . 
and (fT5|) constitute our new extended-MCT equations 
for the coherent density correlator 4)q{t); corresponding 
equations hold for the tagged-particle density correlator 
4>q{t) with the aforementioned replacement of 4>q, rnq, fi^. 



nid 



and 6q by 



i^V 



s\2 



5 id 



and 5q, respectively. 



The extended-MCT equations for the mean-squared dis- 
placement and the non-Gaussian parameter, which are 
required for our discussion in Sec. Illli are derived in Ap- 
pendix [B] All these equations can be solved provided Sq 
and A^ are known as input. (We notice that /„ and /| 
can be obtained based on the knowledge of 5', |10l|.) 



III. RESULTS AND DISCUSSION 

In the following, numerical results of the extended the- 
ory will be presented for the Lennard-Jones (LJ) sys- 
tem in which particles interact via the potential V{r) — 
4eLj{(o'Lj/?')^^ — (ctlj/?")^}. Sq shall be evaluated within 
the Percus-Yevick approximation [22]. This model has 
been studied in Ref. [2^ based on the idealized MCT. 
From here on, all quantities are expressed in reduced 
units with the unit of length ctlj , the unit of energy elj 
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FIG. 1: (Color online) Coherent density correlators (fiq{t) as a 
function of logjo t for reduced temperatures e = {Tc—T)/Tc = 
-0.10, -0.05, -0.03, -0.01, +0.01, +0.05, and +0.10 (from 
left to right). The wave numbers are q = 7.3 (upper panel) 
and q = 10.0 (lower panel), which correspond to the first-peak 
and first-minimum positions of Sq, respectively. The solid 
and dashed curves denote the results from the extended and 
idealized MCT, respectively. The arrows in the upper panel 
refer to the peak positions of the non-Gaussian parameter 
a2{t) for e = +0.01, +0.05, and +0.10 (c/. Sec HITOIl . 



The dynamics as a function of T shall be considered for 
a fixed density p = 1.093, for which the critical tempera- 
ture of the idealized MCT is found to be T^ « 1.637 ^. 
For A^ entering into the extended-MCT equations, we set 
A^ = 0.10 estimated in Ref. !17|| from migration proper- 
ties of crystals. We notice that this value of A^ is consis- 
tent with the Lindemann length [33] ■ We also confirmed 
that the results to be presented below do not rely on the 
specific value of A^ = 0.10: nearly the same results were 
obtained with other values of A^, as far as those values 
consistent with the Lindemann length are chosen. 



A. Coherent density correlators 

Figure [T] shows the coherent density correlators (j)q (t) 
for representative reduced temperatures e = (Tc — T)/Tc 
whose values are specified in the caption. The wave num- 



bers shown are q = 7.3 (upper panel) and q = 10.0 
(lower panel), which correspond to the first-peak and 
first-minimum positions of 5^, respectively. The dashed 
curves refer to the idealized-MCT results which exhibit 
the ergodic (e < 0) to nonergodic (e > 0) transition at 
T = Tc (e == 0) [l3, 111. The solid curves denote the 
results from the extended MCT. 

It is seen from Fig. [1] that the solid curves for e — —0.10 
and —0.05 are hardly affected by the hopping processes, 
but a slight deviation from the dashed curve is dis- 
cernible in the a-relaxation regime of the solid curve for 
e = —0.03. The solid curve for e = —0.01 exhibits the 
same decay up to log^Q i ?» 2 as the corresponding dashed 
curve, but the relaxation thereafter is considerably accel- 
erated. The effects from the hopping processes are drastic 
for e > where the idealized MCT predicts the arrested 
dynamics at the plateau /g, whereas the corresponding 
solid curves from the extended theory relax to zero for 
long times. 

It would be interesting to analyze these extended- 
MCT results based on various scaling laws developed in 
Ref. [ij]. Such an analysis, however, shall be deferred 
to subsequent publications, and in the following, we will 
focus on the connections of the hopping processes with 
the breakdown of the SE relation and with aspects of 
dynamical heterogeneities. 



Breakdown of the Stokes-Einstein relation 
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FIG. 2: (Color online) Logarithmic representation of the 
product Drq* as a function of the reduced temperature e = 
(Tc — T)/Tc. The solid and dashed curves, both referring 
to the left scale, denote the results from the extended and 
idealized MCT, respectively. Circles (right scale) denote the 
Dr]/T data of salol taken from Ref. [l|] and plotted versus e 
with Tc = 262.7 K ^. The arrow marks Tg = 218 K of salol. 
Both the left and right ordinates range over 1.3 decades. The 
inset exhibits the ratio R of the product Drq* to the one at 
a reference temperature where the SE relation holds. The 
prediction from the extended MCT (solid curve) is compared 
with simulation results for a binary mixture of Lennard- Jones 
particles [33| (filled squares) and a binary mixture of dumbbell 
molecules of elongation ^ = 0.8 [33] (filled triangles) . The 
procedure for the comparison is detailed in Ref. [39|. 



Here we investigate the breakdown of the SE relation 
based on the extended MCT, and compare our theoret- 
ical prediction with related computer-simulation results 
and experimental data. A convenient experimental mea- 
sure of the breakdown is a product Drj/T formed with 
the diffusion coefficient D and the viscosity rj, which 
grows as the SE relation fails [i|. In the present study, 
the a-relaxation time r^. of the coherent density corre- 
lator at the peak position q* = 7.3 of Sq, defined via 
the convention 4>q*{Tq*) — 0.1, shall be used as a substi- 
tute for rj/T . This is justified since the T dependence 
of the a-relaxation time at the structure factor peak 
is known to track that of rj/T [3J, [3a]. The diffusion 
coefficient is determined from the long-time asymptote 
D = linit^oo 5r'^{t)/Qt of the mean-squared displacement 
(5r2(i) = {[r,{t) - fs(0)]2), whose extended-MCT equa- 
tions are derived in Appendix IB II 

Figure [2] shows the theoretical prediction for the prod- 
uct DTq* as a function of the reduced temperature e = 
(Tc — T)/Tc. The solid and dashed curves, both referring 
to the left scale, denote the results from the extended and 
idealized MCT, respectively. The idealized-MCT resuff 
for DTq* varies little for e < (the increase is only about 
10% for the e range shown in Fig. [2) , and does not ac- 
count for the breakdown of the SE relation. This reffects 
the universal a-scale coupling predicted by the idealized 
MCT [10] , according to which both the a-relaxation time 
Tq* and the inverse of the diffusivity 1/ D exhibit a uni- 



versal power-law behavior ] e ] '''for small e, and hence, 
the product Drq- approaches a constant for e — > 0— 
(T -> Tc-f). On the other hand, the extended MCT 
predicts the increase of the product DTq* for e ^ —0.05. 
Figure [H thus shows one of the main results of this pa- 
per that the hopping processes are responsible for the 
breakdown of the SE relation. 

Also presented in Fig. [5] is the comparison of the the- 
oretical result with the Drj/T data of salol taken from 
Ref. 1] (circles, right scale). The experimental data are 
also plotted versus e with Tc = 262.7 K of salol deter- 
mined in Ref. [36| . A direct comparison between the the- 
oretical result and the experimental data can be made by 
plotting them on logarithmic scales of the same range as 
done in Fig. [2] It is seen that the extended-MCT result 
is consistent with the experimental data concerning the 
degree of the breakdown of the SE relation. 

The inset of Fig. [2] compares the extended-MCT re- 
sult (solid curve) with simulation results for a binary 
mixture of Lennard- Jones particles [37[ (filled squares) 
and a binary mixture of dumbbell molecules of elonga- 
tion C = 0.8 [38] (filled triangles). Here the comparison 
is done in terms of the ratio R of the product DTq<. to 
the one at a reference temperature where the SE relation 
holds. (See Ref. [3^ concerning the details of the com- 
parison.) By definition, the ratio R is unity if the SE 
relation holds, whereas it exceeds unity as the SE rela- 



tion fails. Again, the theoretical prediction is consistent 
with the simulation results. 

On the other hand, recent measurements [2,y| indicate 
that the self-difFusion coefRcient is about 100 times faster 
near Tg than that predicted by the SE relation. This is 
about a factor of 10 larger compared to our theoretical 
prediction near Tg {cf. Fig. [2]), assuming that Tg of the 
present system is located at e « 0.17 estimated from the 
reduced temperature at Tg of salol. This implies that our 
model for the hopping kernel might be too primitive to 
be applicable near Tg. In the following, we shall therefore 
focus mainly on the regime e ^ +0.05 where our theoret- 
ical prediction is consistent with the simulation results 
and experimental data. 




C. Aspects of dynamical heterogeneities 

1. Non-Gaussian parameter 

We next explore a connection of the hopping pro- 
cesses with aspects of dynamical heterogeneities. We 
start from the discussion on the non-Gaussian parameter 
CK2(t) which characterizes deviations from the Gaussian 
behavior of the van Hove self correlation function [401 : 



a2it)^-[Sr^it)/Sr^{tr]-l. 



(16) 



Here 6r'^{t) = {[rsit) — f'siO)]'^)- General properties of 
a2{t) revealed b y c omputer simulations can be summa- 
rized as follows [l3|: (i) on the time scale at which the 
motion of particles is ballistic, a2{t) is zero; (ii) upon en- 
tering the time scale of the /3 relaxation where the density 
correlators are close to their plateaus, a2{t) starts to in- 
crease; and (iii) on the time scale of the a relaxation, 
a2{t) decreases to zero, reflecting diffusive dynamics at 
long times which is a Gaussian process. It is also observed 
that the maximum value of a2 (t) and the time at which 
this maximum is attained both increase with decreasing 
T [1^. Positive a2it) means that the probability for a 
particle to move very far is enhanced relative to the one 
expected for a random-walk process. The peak height 
of a2{t) has therefore been interpreted as a measure of 
the dynamical heterogeneity that reflects different local 
environments around an individual particle. 

The extended-MCT equations for determining the non- 
Gaussian parameter are derived in Appendix lB 2\ and the 
resulting a2 (t) are plotted in Fig. [3] for reduced temper- 
atures e = (Tc — T)/Tc specified in the caption. The 
corresponding results from the idealized MCT, but for 
e < (i.e., T > Tc) only, are presented as dashed curves 
in the inset to avoid the overcrowding of the figure. As 
noticed at the end of Appendix IB 21 the peculiar behav- 
ior a2{t — > 0) = —2/3 for short times as predicted by 
our theory simply reflects that the ideal-gas contribution 
to the memory kernel is discarded in the mode-coupling 
approach. Such an ideal-gas contribution is responsible 



FIG. 3: (Color online) Non-Gaussian parameter 0:2 (i) as a 
function of logj^Q t based on the extended MCT for reduced 
temperatures e = (Tc - T)/Tc = -0.10, -0.05, -0.03, -0.01, 
-1-0.01, -1-0.03, and +0.05 (from left to right). The correspond- 
ing results from the idealized MCT, but for e < (i.e., T > Tc) 
only, are presented as dashed curves in the inset to avoid the 
overcrowding of the figure. Notice about an order of magni- 
tude difference in the ordinate scales. 



for the short-time ballistic regime, but is irrelevant as far 
as the long-time dynamics is concerned. 

It is seen from the inset of Fig. [3] that the peak height 
of a2{t) predicted by the idealized MCT does not grow 
with decreasing T, and it is underestimated by almost an 
order of magnitude compared to that reported in com- 
puter simulations [19|. This defect has already been 
known from the theoretical work in Ref. [S ill, and has 
been demonstrated more explicitly in Ref. [2l| . The main 
panel of Fig. [3l on the other hand, indicates that a sub- 
stantial improvement is achieved by the extended MCT 
in that the peak height grows upon lowering T to an ex- 
tent as observed in simulations [10|. (Notice an order of 
magnitude difference in the ordinate scales of the main 
panel and the inset.) Thus, the extended MCT repro- 
duces a feature of the dynamical heterogeneity charac- 
terized by the non-Gaussian parameter, and this is ac- 
complished via the inclusion of the hopping processes. 



2. Probability distribution of particle displacements 

In terms of the non-Gaussian parameter a2(i), the dy- 
namics is most heterogeneous in the late-/? or early a 
regime where the density correlators start to decay from 
the plateau (see the upper panel of Fig. [1]). Recently, 
it has been recognized from studies of four-point density 
correlation functions that there exists dynamical hetero- 
geneity on a much lon ger time scale comparable to the 
a-relaxation time r^* [4l|. Recognizing that a2{t) is 
dominated by those particles which move farther than 
expected from a Gaussian distribution of particle dis- 
placements, Flenner and Szamel introduced a new non- 
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FIG. 4: Probability distribution P(logjQ (5r);t) of the loga- 
rithm of single particle displacements based on the extended 
MCT at the reduced temperature e = (Tc - T)lT^ = +0.01 
(solid curves). The times shown are f = 5.4 x 10^, 4.3 x 10"', 
1.7 X 10^ 3.4 X 10^ 6.8 x 10^ and 2.7 x 10^ (from left to 
right). The Q-relaxation time r,* at this reduced tempera- 
ture is Tq» — 4.3 x 10^ (c/. Fig. [l]). The dashed curves refer 
to the probability distribution under the diffusion approxi- 
mation for i = 6.8 X 10^ and 2.7 x 10®. The inset exhibits 
the time evolution of the mean-squared displacement 5r'^(t) 
at e = -1-0.01 on double logarithmic scales (solid curve). The 
dotted curve refers to the diffusion asymptote, 6Dt. 



Gaussian parameter 7(i) = (l/3)(5r^(t)[5r^]~-'^(i) — 1 with 
\8r'^\~^[i) = (l/[rs(i) — rs(0)]^), which weights strongly 
the particles which have not moved as far as expected 
from the Gaussian distribution 20]. It is found that the 
peak position of 7(i) is located on the time scale of r^*, 
and hence, characterizes the longer-time dynamical het- 
erogeneity. It is also observed that the peak position 
of 7(t) corresponds to the time at which two peaks in 
the probability distribution of particle displacements, re- 
flecting populations of mobile and immobile particles (see 
below), are of equal height. This implies that the pres- 
ence of the longer-time dynamical heterogeneity can be 
examined also through such a probability distribution. 

The mentioned probability distribution P(logiQ(5r); t) 
of the logarithm of particle displacements 6r at time t can 
be obtained from the van Hove self correlation function 
Gs ((5?", t) by the transformation [23| 

P(logio(fr); t) = In(lO) 4^fr3G,(fr, t). (17) 

The probability distribution is defined such that the inte- 
gral / ^ dx P(x\ t) is the fraction of particles whose value 
of log]^g((5r) is between xq and xi. If the motion of a 
particle is diffusive with a diffusion coefficient I), there 
holds Gs{5r,t) « [l/{4TTDtf^^]exp{-6r^/ADt) ^. As 
argued in Ref . [2d| , the shape of the corresponding proba- 
bility distribution under the diffusion approximation be- 
comes time independent with the peak height w 2.13. 

Figure[3]shows the extended-MCT result for the proba- 
bihty distribution P(\ogiQ{6r);t) (solid curves) at the re- 
duced temperature e = (Tc - T)/Tc = +0.01. The times 



shown are i = 5.4 x 10^, 4.3 x 10^, 1.7 x 10^, 3.4 x 10^ 
6.8 X 10^, and 2.7 x 10^. The a-relaxation time r,. at 
this reduced temperature is r,. — 4.3 x 10^ (c/. Fig.[T]). 
The dashed curves refer to the probability distribution 
under the diffusion approximation for t — 6.8 x 10'^ and 
2.7 X 10^. The inset exhibits the mean-squared displace- 
ment dr'^{t) at e = +0.01 on double logarithmic scales, 
from which one understands that the times chosen in the 
main panel range from the late-/3, plateau regime to the 
final a regime where 6r'^{t) ~ 6Dt. 

The appearance of the plateau in Sr^ (t) is due to par- 
ticles being caged [31|, and the peak of the probability 
distribution P{\ogiQ{6r);t) for i = 5.4 x 10'^ reflects pop- 
ulations of such "immobile" particles. At later times, 
the second peak develops in the probability distribution 
at \ogiQ{Sr) « 0.0 (i.e., 6r « 1.0), reflecting "mobile" 
particles hopping over interparticle distances. One in- 
fers from Fig. [5] that the time at which the two peaks in 
P {\ogiQ {Sr);t) become of equal height is located on the 
time scale of r^* , and this is consistent with the simula- 
tion result [201 • Subsequently, the double-peak structure 
disappears, and the probability distribution approaches 
the one well described by the diffusion approximation. 

Coexistence of mobile and immobile particles is a di- 
rect indication of the dynamical heterogeneity [4i]. Our 
results shown in Fig. [5] indicate that the hopping pro- 
cesses are responsible for such a double-peak structure in 
the probability distribution P{\og^Q{Sr);t) occurring on 
the time scale Tg* of the a relaxation. This also explains 
why the idealized MCT failed to rep roduce the double- 
peak structure in P{\ogiQ{Sr);t) J2l| . 



3. Growing dynamic length scale 

Figure S] also implies that the probability distribution 
approaches its diffusion asymptote only after its peak po- 
sition exceeds those length scales where the double-peak 
structure in P{\ogiQ{6r);t) is observable. Thus, there is 
a certain length scale above which Fickian diffusion sets 
in. In the following, we shall quantify such a length scale 
characterizing the crossover from non-Fickian to Fickian 
diffusion, and investigate its temperature dependence. 

To this end, let us introduce the ratio R^ of the prod- 
uct q^Dr^ to the one at some reference temperature, 
which is an analogue of the ratio R studied in the inset 
of Fig. [2l Here the a-relaxation time r^ of the tagged- 
particle density correlator is deflned via the convention 
(f>i{T^) = 0.1. Under the diffusion approximation, there 
holds (?!)^(t) w exp{-q'^Dt) ^. Thus, the product q'^Dr^ 
is constant, and hence, the ratio R^ is unity, if the dynam- 
ics on the length scale « ^ir/q is diffusive. On the other 
hand, the ratio R^ exceeds unity if the dynamics on the 
length scale w ^ir/q is non-Fickian. We shall therefore 



viai?^ 



1.1, 



define the crossover wave number gonsct 
i.e., as the wave number at which the ratio R^ reaches 
10% above unity. The onset length scale of Fickian dif 
fusion shall then be defined via ^onsct = 27r/(7onsct • 
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FIG. 5: (Color online) Ratio Rq of the product q^Dr^ to 
the one at a reference temperature Tref as a function of the 
reduced temperature e = (Tc — T)/Tc. Here Trcf is chosen 
such that (Tc — T,.ai)/Tc — —0.125. The wave numbers shown 
are q = 7.3 (solid curve), 4.0 (dashed curve), 2.0 (dashed- 
dotted curve), and 1.0 (dotted curve). The horizontal dashed 
line marks Rg = 1.1 chosen to determine the crossover wave 
number gonsot introduced in the text. The inset exhibits the 
onset length scale of Fickian diffusion ^onset = 2n/qonset as a 
function of the reduced temperature e. 



Figureiniexhibits the ratio R^ for wave numbers q — 7.3 
(solid curve), 4.0 (dashed curve), 2.0 (dashed-dotted 
curve), and 1.0 (dotted curve), as a function of the re- 
duced temperature e = (Tc — T)/Tc. The reference 
temperature Trcf for calculating R^ is chosen such that 
(Tc - T,^{)/Tc = -0.125. The horizontal dashed line 
marks R^ — 1.1 chosen to determine the crossover wave 
number Qonsct introduced above. 



For reduced temperatures e < —0.04, all the ratios 
R^ shown in Fig. [5] are less than 1.1, implying that the 
tagged-particle dynamics on the length scales comparable 
to and larger than 27r/7.3 w 0.86 can be well described by 
the Fickian diffusion law. For e > —0.04, R^ for q — 7.3 
exceeds 1.1, meaning that the dynamics is diffusive only 
on length scales larger than 0.86. The onset length scale 
^onsct of Fickian diffusion increases to 27r/4.0 « 1.6 at 
e « —0.03, and to 27r/2.0 « 3.1 at e « 0, as can be 
inferred from Fig. [51 The length scale ^onsot so obtained 
as a function of the reduced temperature e is summarized 
in the inset of Fig. [5] Thus, the extended MCT predicts 
the presence of a growing dynamic length scale. It is 
anticipated that ^onsct is intimately related to the mean 
size of the dynamic clusters observed in simulations [g] , 
since Fickian diffusion, i.e., a random- walk process, is 
possible only over those length scales where the coherence 
length of such clusters is smeared out. 



IV. SUMMARY 

In this paper, we developed a new extended version of 
MCT for glass transition. The activated hopping pro- 
cesses are incorporated via the dynamical theory, origi- 
nally formulated to describe diffusion-jump processes in 
crystals and adapted in this work to glass-forming liquids. 
The dynamical-theory approach treats hopping as aris- 
ing from vibrational fluctuations in quasi-arrested state 
where particles are trapped inside their cages, and the 
hopping rate is formulated in terms of the Debye- Waller 
factors characterizing the structure of the quasi-arrested 
state. The resulting expression for the hopping rate takes 
an activated form, and the barrier height for the hopping 
is "self-generated" in the sense that it is present only in 
those states where the dynamics exhibits a well defined 
plateau. It is discussed how such a hopping rate can be 
incorporated to develop a new extended MCT. 

The extended MCT deals with the interplay of two ef- 
fects. Nonlinear interactions of density fluctuations, as 
described by the idealized memory kernel m^^{z), lead to 
the cage effect with a trend to produce arrested states for 
sufficiently low temperatures. Phonon assisted hoppings, 
taken into account via the hopping kernel Sq(z), lead to 
the a relaxation at all temperatures and restore ergodic- 
ity. The interplay of these two effects is described by the 
memory kernel mg{z) via Eq. (j3bp which has the form 
of a Dyson equation. As demonstrated in Sec. IIIII for a 
model of the Lennard- Jones system, it leads to nontriv- 
ial theoretical predictions concerning the breakdown of 
the SE relation and characteristic features of dynamical 
heterogeneities, which the idealized version of the theory 
failed to reproduce. 

In dense liquids, relaxation is necessarily connected 
with rearrangements of large complexes of particles. We 
have seen in Fig. [3] that the peak height of the non- 
Gaussian parameter a2 (t) grows substantially in the late- 
f3 regime as the temperature is decreased. This enhanced 
probability for a particle to move further is what one 
would expect as a result of the building of a backflow 
in the liquid, and it is anticipated that the string-like 
motions observed in the late-/3 regime [i^l reflect such 
a backflow pattern. The backflow was originally dis- 
cussed by Feynman and Cohen and taken into account 
in their theory for roton spectrum in liquid helium [43l |. 
Subsequently, it was found that the quantum-mechanical 
analogue of the idealized memory kernel re pro duces the 
Feynman-Cohen result for roton spectrum [44| . i.e., the 
backfl ow p henomenon is within the scope of the idealized 
MCT [2^. However, the idealized memory kernel alone 
does not account for such a pronounced peak in a2{t) as 
observed in computer simulations [315. This implies that 
the mentioned interplay with the hopping kernel plays a 
relevant role for the substantial increase of a2 (i) and the 
building of the backflow. 

The double-peak structure in the probability distribu- 
tion P{\ogiQ{dr);t) of particle displacements, which is 
most pronounced on the time scale t„* of the a-relaxation 



and disappears at longer times (c/. Fig. |4]), reflects co- 
existence of mobile and immobile particles with a life 
time « Tq- , and is a direct indication of the dynamical 
heterogeneity [4|. The extended MCT developed here 
provides a natural explanation for its origin in terms of 
the cage and hopping effects. As argued in Sec. IIII C 3) 
the life time of the double-peak structure gives rise to a 
growing dynamic length scale £onsot, which is associated 
with the decoupling of the time scales of the a processes 
occurring on length scales smaller than ^onsot from the 
diffusion coefhcient. The breakdown of the SE relation 
discussed in Sec. IIIIBI can also be understood in terms 
of the decoupling of the a-relaxation time Tq* from the 
diffusivity which occurs when ^onsct exceeds the average 
interparticle distance « 27r/g*. It is anticipated that the 
onset length scale fonsct of Fickian diffusion is intimately 
related to the mean size of the dynamic clusters observed 
in simulations [6i], since Fickian diffusion, i.e., a random- 
walk process, is possible only over those length scales 
where the coherence length of such clusters is smeared 
out. Then, our picture for the breakdown of the SE re- 
lation is consistent with that of Ref. [0| , where the con- 
nection between decoupling phenomena and the growing 
coherence length scale is discussed. 

It is certainly necessary to improve further the theory 
developed here. First of all, possible correlated hopping 
effects are discarded in our model for the hopping ker- 
nel. This might explain why our theory underestimates 
the degree of the breakdown of the SE relation near Tg, 
which we referred to in connection with Fig. [2l That the 
development of the length scale ^onsot seems suppressed 
for e > 0, as can be inferred from the inset of Fig. [5l 
might also be related to such a defect of the present the- 
ory. Second, we did not determine A^ - a square of the 
ratio A = x* /a formed with the critical size x* of the 
phonon-assisted fluctuation needed to cause a hopping 
and the average interparticle distance a {cf. Sec. |ll]) ~ 
microscopically, but simply took a value (A^ = 0.10) 
from the literature. Regarding this, let us mention that 
nearly the same results as those presented in Sec. lIIII can 
be obtained with other values of A^ , as far as those val- 
ues consistent with the Lindemann length 32] are chosen. 
Thus, our theoretical results do not rely on the specific 
value of A^ = 0.10. But, of course, it is desirable to 
determine A^ consistently within the theory. Third, the 
a-relaxation stretching is not enhanced at low temper- 
atures. For example, the stretching exponents f3q, ob- 
tained via Kohlrausch-law fits of the density correlators 
shown in Fig. [H for q = 10.0 are 0.76, 0.79, and 0.80 for 
e = -1-0.01, -1-0.05, and -t-0.10, respectively. One possi- 
bility to overcome this defect is to take into account dis- 
tribution of hopping rates or barrier heights, which can 
be done, e.g., by considering fluctuations in the Debye- 
Waller factors and A^. Four-point density correlators 
might be necessary for this purpose, though it is a diffi- 
cult task to calculate such higher-order correlators. But, 
in view of the significant results achieved by our theory, 
it is promising to pursue its further development. 
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APPENDIX A: NOVEL DERIVATION OF THE 
HOPPING KERNEL OF GOTZE AND SJOGREN 

In order to go beyond the idealized MCT, one has 
to consider corrections to the idealized memory kernel, 
mq{z) = m^(z) + Amq{z). The extended MCT of Gotze 
and Sjogren [Ij, |23| can be considered as a theory for 
such corrections, and yields an expression for Amq{z) 
in terms of the hopping kernel Sq{z) that involves cou- 
plings to currents. Their extended theory is formulated 
with the generalized kinetic theory for phase-space den- 
sity fluctuations, but essentially the same expression for 
the hopping kernel Sq{z) can be derived also with the 
standard projection-operator approach based on density 
and current-density fluctuations. In this appendix, we 
outline such a novel derivation. 

We start from reviewing the derivation of the ideal- 
ized memory kernel mjf (t) (see Ref. [Ifll for details). Let 
us introduce a projection operator V onto the subspace 
spanned by the density fluctuations p^ and the current- 
density fluctuations j^ = ^,- v^ exp(z(f ■ r^) (A — x, y, z). 

Here vf denotes the A component of the velocity of «th 
particle. Within the Zwanzig-Mori formalism [23], one 
obtains based on the operator V the exact equation ((Ta|) 
for the coherent density correlator 4>q{t). The formal ex- 
pression for the memory kernel entering there reads 



nlruqit) 



M 



NksT 



(^1 



JQCQt 



Fq), 



(Al) 



where Q = 1 — V, and the fluctuating force is given by 

(A2) 



Fq = dtiq-j^) - iq--r^Pq. 



'MSa 



Here q — q/q. Under the mode-coupling approach, the 
fluctuating forces are approximated by their projections 
onto the subspace spanned by pair-density modes /c^rPp; 
this is done by introducing the second projection operator 



v.x 



1 



2^PkPp^2s,s, 



(Plpp^)^ 



(A3) 



k>p 



and approximating F^ « 7^2 -^V = F^'' ^ i.e.. 



"«™'('^ "" ]^(^r*e'^''''*^f )• (A4) 
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Notice here that the static factorization approximation 

iPlPfPk'Pp') « h u'^P.P'N^SkSp {k > p, k' > p'), (A5) 



is aheady used in the definition (|A3|) of V2 ■ Within the 
convolution approximation for triple density correlations, 
one finds for the projected fiuctuating force 



F^ 



pp 



. pk-gT 

'* NM 



/ ' q,k+p 

k>p 

x{q-kck + q- pcp)pj:pj; (A6a) 

(A6b) 



= -'l^]^2^1-^'"'PkPg-k 



Substituting Eq. (|A6aP into Eq. (jA4p and then using the 
dynamical factorization approximation 

iplPpe^'^^'^'ppPp') « Sf:^pS,~p^,N'S,SpMt)Mt) 

{k>p,k'>p'), (A7) 

which factorizes averages of products evolving in time 
with the generator Q£Q into products of averages formed 
with variables evolving with C, one obtains the idealized 
kernel m]^{t) given in Eqs. ^. 

To extend the idealized theory, one has to avoid the 
use of the dynamical factorization approximation (|A7[) . 
We shall therefore start from Eqs. (|A4p and (IA6bp for the 
memory kernel mg{t), to which the approximation (|A7[) 
has yet to be applied. Notice, however, that the static 
factorization approximation (jASp is still assumed in our 
approach; this is employed in the definition (|A3[) of 7^2, 
which is then used in deriving Eq. (|A6b[l . Thus, there 
holds mg(0) = mf{0) at i = 0. 

By introducing a new projection operator V' onto the 
subspace spanned by F£'^, i.e., 

V'X = F£P ^ /FPP*X), (A8) 

one obtains the following Zwanzig-Mori equation of mo- 
tion for mq{t): 

dtmq{t) + m,(0) / dt' Lq{t - t')mg{t') = 0. (A9) 



Here the factor mq{0) in front of the convolution integral 
is just for later convenience, and the formal expression 
for the "memory kernel" Lq(t) is given by 



mq{0)Lq{t) 



M 



{R}e 



NkBTQ^qlTlqiO) 

in terms of the "fluctuating force" 



iQ'QCQQ'tj^^^^ 



(AlO) 

(All) 



where we have introduced Q' = 1 — V' . Substituting 
Eq. ([Mb)) into Eq. (|MT|) . one finds 
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pksT 
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k 

-{k-Jq^Sp+p-Jq-Sk)}, (A12) 



in which p = q — k. 

Now, we apply the mode-coupling approximation to 
the kernel Lq{t). Because of the odd time reversal sym- 
metry of R^, the simplest mode-coupling approximation 
for Lq{t) can be introduced by defining the projection 
operator V2 onto the subspace panned by the product 
PjTji of the density and current-density modes: 



nx = Y^p,3^ 



M 



N-^k^TSk 



{plji*X), (A13) 



where we have used the factorization approximation 

iP*Pf*Pk'Jfi') « Sj:j:,6^^N^kBTSk/M. (A14) 

Here and in the following, we use abbreviations p = q — k 
andp' = q — k'. We thus obtain under the mode-coupling 
approximation R^ w ^2-^? = R'^' 

mq{0)Lq{t) 



<-q ~ r 2-'<-(j q 

M 



lrjP3*ptQ'QCQQ'tTjpjx 

NkBTnlmqiOy^'T ^^'■ 

(A15) 
Within the convolution approximation for triple density 
correlations, one finds 

k ^ 



in which V\ is given by 



V\{q;k,p) = 



dk' q ■ k' Ck' 



(27r)3 
X (fcl^p'/i|p^-_£| +p'xSk'h^i:,_j:^). (A16b) 

Here kx refers to the A component of the vector k, and the 
function hq is given by hq — CqSq. When Eq. (|A16a|) is 
substituted into Eq. (|A15|) . the kernel Lq{t) is expressed 
in terms of four-mode correlators, for which we invoke 
the following dynamical factorization approximation: 

^h,M''''fiW^'^'Pv) 
^6^pN^ikBT/M)SkMt) 



6j:^^,N^SkSpCk^Px/kp)dtMt)dtMt)-i^^7) 
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Here </>Wi) and (/)J{t) are longitudinal and transversal 
current correlators, respectively, which are normalized to 
unity at i = 0. We then obtain the following form for the 
kernel Lq{t) that involves couplings to current modes: 



+ VT{q;k,p)Mt)(f^it) 

+ V'{q- k,p)dtMt) dtMt)] ■ (Al8a) 

The vertex functions Vl, Vt, and V are expressed in 
terms of the thermal velocity, the static equilibrium 
quantities, and Vx given in Eq. (jA16bp as 



VL{q;k,p) = 



VT{q;k,p] 



pksT SgS, 



q'Jk 



{2tt)^M q^TUqiO)^ 



^^PxVxiq;k,p) , 

(A18b) 



-{Y,V^i<l:k,pf 



pksT SqSk 
{2nfM q^mq{0)^ 

-[Y,PxVx{q;k,p)f}, (Al8c) 



V'{q;k,p) = 



t>qJk^' 



(2^)3(?2^,(0) 



f,[ET^A(g;p,fc) 



\Y.^vx{q-Xp: 



(A18d) 



Finally, let us connect the Laplace transform of Lq {t) 
to the hopping kernel Sq{z) using the definition given by 
Eq. ([5a|) . The Laplace transform of Eq. (|A9[) reads 



rriqiz) = - 



rUqiO) 



z + mq{0)Lq{z) 



(A19) 



In view of this expression, let us formally introduce the 
function L'f^{z) in terms of the Laplace transform of the 
idealized memory kernel m^^(z) via 



m;<^(z) 



m,(0) 



z + mq{0)Dq^iz)' 



(A20) 



Here we have used the equality mq{0) = ttt,' (0) at < = 
noticed above. Substituting Eqs. (|A19p and (|A20p into 
Eq. (I3ap . one obtains for the hopping kernel 



Sqiz)=Lq{z)-Vq-iz) 



(A21) 



From the functional form of Lq{t) given in Eq. (jA18a[) 
and interpreting that L'-^{z) in Eq. (|A2ip subtracts 
those contributions already accounted for by the ideal- 
ized memory kernel m^^ (z) , one understands that the ex- 
pression (|A2ip for the hopping kernel is essentially the 
same as the one derived by Gotze and Sjogren |13|.|23|. 



APPENDIX B: EXTENDED MCT EQUATIONS 

FOR THE MEAN-SQUARED DISPLACEMENT 

AND THE NON-GAUSSIAN PARAMETER 

1. Mean-squared displacement 

The equation of motion for the mean-squared displace- 
ment Sr'^(t) can be obtained from Eq. (|lap for 4>q{t) by 
exploiting its relation to the small-^/ behavior of 0* (t) = 
l~q^6r\t)/6 + Oiq'^) \^: 



dt5r\t)- 



k^T 

M 



dt'mlg){t-t')dr^{t') 



,'-§t. m 



Here we have introduced the g — > limit of the memory 
kernel via ml^Jt) = limg_,o q^mq{t). To obtain ra'^,^At), 
it is more convenient to rewrite Eq. (|3bp in the form 



m^(z) = mf\z) + m^'^(z) <5^(z) m^(z), (B2) 

from which one finds 

^Ui^) = ^wi-A + ^l^i^) 5(0)(^) "^(0)(^)- (B3a) 

Here m?Q?(z) refers to the corresponding q ^ limit of 
the idealized memory kernel, and (5?„s (z) is defined via 
Jf^Jz) = hmq^o 5q{z)/q'^- The former is given by [3i] 

™(o?(0 = i / ^fc k^pSkclc^k{tWk{t), (B3b) 

whereas the latter reads after carrying out the g ^ 
limit in Eq. dH]) 



^(o)(^) = iwiiopNcO^ /Q. 



(B3c) 



Equations pjlj) and p33)) constitute the extended-MCT 
equations for the mean-squared displacement 6r^ {t) . 



2. Non-Gaussian parameter 

The non-Gaussian parameter a2{t) defined in Eq. (fT6|l 
can be obtained from the mean-squared displacement 
5r^{t) and the mean-quartic displacement 5r'^{t). The 
extended-MCT equations for 5r^ (t) can be derived using 
the same method employed above for Sr'^{t), but with 
higher order expansions in q. 

Since Sr'^{t) is proportional to the fourth Taylor coef- 
ficient in the small-g expansion of 4>q{t) = 1 — q^6r^ /3l + 
q'^Sr'^/bl + 0{q^) Q, one can derive the following equa- 
tion from the small-q behavior of Eq. (|lap for 0^ [t] : 



dtdr^{t)-20{kBT/M) dt' 6r^{t) 
Jo 
ft 



+ {kBT/M) dt' m-L.{t-t')Sr^{t') 
1 L 

10Tnl2){t-t')Sr^{t')] =0. (B4) 
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Here we introduced the memory kernel mf„> (t) via the find for short times 



small-g expansion of q^m^{t) 



'(0) 



it) 






(2) 



(t)/2. 



0(g'^). Using the corresponding expansion for the 
ideahzed memory kernel m'l^'^{t) and the expansion 
5l{z)/q^ = Sf^^iz) + qHl^^{z)/2 + 0{q^) for the hopping 



kernel, one obtains from Eq. (|B2p 

+ "^(2)W'5^o)W"^(o)W- (B58 
The expression for m'f^) (t) reads [3i| 
1 



'(2) 



(t) 



107r2 



(9fc2 3fc 9A: 



, (B5b) 



while one obtains from the small-g expansion of Eq. ([T 



(5(^2) (^) ^iw^opN,a'[ri ^ a' /2Q\/i 



(B5c) 



Here, r^ is defined via the small-q expansion of the Lamb- 
Mossbauer factor /| = 1 - q'^rl + 0(g'') [si^. Equa- 
tions (JB4|) and (JBSJl along with Eqs. (JB3|l constitute the 
extended-MCT equations for (5r^(i). The non-Gaussian 
parameter a2{t) can then be obtained from Eq. pB|) . 

Let us consider the short-time behavior of a2 [t) based 
on the derived equations. From Eqs. (|B1|) and (jB4p we 



5r2(i) = 3(fcBr/Af)2t2 + 0(i4)^ 



1 



'(2) 



(0) 



(B6) 

i4 + 0(f6).(B7) 



5r^{t) = b{kBT/Mf 
Substituting these results into Eq. (I16p . one obtains 

-"^(2)(0) 



a2(t) = 



6 



0{i' 



(B8) 



Since 0^(0) 



1, we obtain from Eq. (|B5b)) mUK'i) = 



Using this result, one can show that m%AQ) = based 
on Eqs. (|B5ap and (|B5cp . According to Eq. (|B8|) . this 
means that the initial value of a2 (t) within the extended 
MCT (and also within the idealized MCT which is based 
on Newtonian dynamics) is given by 



a2{t ^ 0) = -2/3. 



(B9) 



This is in disagreement with the exact initial behavior 
"2 (i — >• 0) = [22 • This discrepancy simply reflects that 
the ideal-gas contribution to the memory kernel is dis- 
carded in the mode-coupling approach. Indeed, the ideal- 
gas contribution yields m'l^JO) = 4 |22l |. which when sub- 
stituted into Eq. jBSj recovers a2(i — ^ 0) = 0. The ideal- 
gas contribution to the memory kernel is responsible for 
the short-time ballistic regime, but is irrelevant as far as 
the long-time dynamics is concerned. 
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